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Far-infrared loops in the 2 nd Galactic Quadrant 
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Abstract. We present the results of an investigation of the large scale structure of the diffuse interstellar medium in the 
2 nd Galactic Quadrant (90° < 1< 180°). 145 loops were identified on IRAS based far-infrared maps. Our catalogue lists their 
basic physical properties. The distribution clearly suggests that there is an efficient process that can generate loop-like features 
at high galactic latitudes. Distances are provided for 30 loops. We also give an observational estimate on the volume filling 
factor of the hot gas in the Local Arm, 4.6% < f 2 mt < 6.4%. 
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1. Introduction 

A study of the diffuse interstellar medium is presented, a step 
in the way deriving the near-past and near-future of the galac- 
tic star formation. The large scale structure of the ISM in the 
Galaxy is diverse with the complex distribution of cavities, fil- 
aments, arc, loops and shell s, which is often re ferred to as the 
"Cosmic Bubble Bath" jBrand & Zeaievlll975l) . Ridges of en- 
hanced radio contiunuum emission extending to high galac- 
tic latitudes have been known since the 60s. The f irst exten- 
■ sively studied examples are Loop ! jLarge et allll966(). Loop II 
("Cet us Arc" Large et all 1 1962b and Loop IIl jLarge et all 
119621) . Later anomalous HI features like worms, chimneys, su- 
pershells were identified and catalogued, up to a character- 
istic size of ~2kpc with expansion velocities of ~20kms~ 1 
dHeilesL 1 19791 1 198(1 1 1984 lllul Il98lt iKoo & Heilesl Il99lh . 
These were as sociated with voids in the HI distribution (see e.g. 
|Palouslfl996l for a review). Similar complex morphology has 
been reco gnized in som e nearby spiral and irregular gala xies 
e.g. M31 (Bri nkslll98ll. M33 JPeul & den HartogMl990l) and 
the M agellanic Clouds istavelev-Smifh et allll992UKim et ail 
1998). 

The most studied galactic shells were formed by su- 
pernova (SN) explosions and win ds of massive stars 
jTenorio-Tagle & Bodenheimeit Il988i). However, high veloc- 
ity clouds (HVCs) may also form large cavities when infalling 
from the galactic halo and colliding with the ISM of the disc. 
A remarkable example of a HVC - galactic disc interaction i s 
the North Celestial Pole (N CP) loop jMeverdierks et allfl99ll) . 
lEhlerova & Paloul i 19961) investigated the origin of galactic 
HI shells and found, that those are likely related to star for- 
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mation, and not to the infall of HVCs. Apart from the con- 
nection of these features to young massive stars and their as- 
sociations, they may also provide a nest for the next genera- 
tion of stars via the process of pr opagating s tar formation . The 
idea was first proposed bv lOortl i 19541) and lOpikl Jl953|). and 
later followed by observational evidences (see Blaauw, 1991, 
for a review). Th i s m echanism was modelled in detail by 
Elmegreen & Lada i 19771) and successfully applied to the se- 
quential formation of subgroups in OB associations. Gamma 
ray bursts may also be responsibl e for the form a tion o f the 
largest HI holes, as proposed by lEfremov et all (1998) and 
lLoeb & PernalJl998l) . Nonlinear development of instabilities in 
the ISM may form large cavities as well, without stellar energy 
injection Jwada et all l2000t iKorpi et all 1 19991 iKlessen et all 
l200rilSanchez-Salcedd.l200ll) . 

The cavities of the I SM contain low density h ot gas in 
the multiphase model bv lMcKee&Ostrikeil < 19771) . The fill- 
ing factor / of this hot gas in the galactic disc is an impor- 
tant param eter of the ISM , whi ch is poorly constrained obser- 
vationallv. iFerrierd i 19981) and iGazol-Patifio & Passoil ^999) 
estimate / m 20% at the solar circle, which drops to / = 3-8% 
for 8.5 <R@< lOkpc galactocentric distances. 

The 21 cm emission of diffuse HI clouds is well-correlate d 
with the 100 ^m FIR emission jBoulanger & PeraultL IT988). 
This exte nded emission is often referred to as the galac- 
tic cirrus JLow et all [T984). A typical far-infrared colour of 
A/60 /AZi oo = 0.25 was derived from IRAS/ISS A maps whe n 
the gains are corrected according to IWheelock et all {T994). 
According to COBE/DIRBE results its dust colour temperature 
is T^ ~ 1 8 K llLagache et allll 998l) . 

Studies of the large scale structure of the cirrus emis- 
sion and loop- or arc-like features were restricted so far to 
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low galactic l atitudes in the far-infrared dSchwartzl Il987t 
llvlarstonl Il996l) . and reached only medium galactic latitudes 
in radio surveys. Detailed investigations of individual loops 
were perfor med only in th e vicinity of the galactic plane. 
Recently Dai gle et"al] J2003) presented the first automatic tool 
for detecting expanding HI shells using artificial neural net- 
works. Their data is restricted to the sky area 75° < 1 < 145°and 
-3° < b < 5°and their detection strategy is based on the velocity 
information of the HI data rather than morphology. 

In order to get a complete, full-scale view on the large-scale 
structure of the cold diffuse ISM, in our study we searched 
the whole 2 nd galactic quadrant (90° < 1 < 180°; -90° < b < 90°) 
for intensity enhancements in the interstellar medium, showing 
loop-like structures. 

2. Data analysis 

2.1. Quest for loops on far-infrared maps 

We investigat ed the 60 and lOO^m ISSA plates (IRAS Sky 
Survey Atlas, IWheelock et all [1.994) in order to explore the 
distribution of dust emission in the 2 nd Galactic Quadrant. We 
created composite images of the 12?5xl2?5 sized individual 
ISSA plates using the "geom" and "mosaic" procedures of 
the IPAC-Skyview package, both at 60 and 100 /mi. These im- 
ages were built up typically from ~ 10-15 ISSA plates, reaching 
a size of ~40°x40°. Loop-like intensity enhancements were 
searched by eye on the 100 /mi mosaic maps. Loops by our defi- 
nition must show an excess FIR intensity confined to an arc-like 
feature, at least 60% of a complete ellipse-shaped ring. A loop 
may consist of a set of bright, more or less isolated, extended 
spots, or may be a diffuse ring or part of a ring. The size of the 
mosaic image limits the maximal size of the objects found. On 
the other hand, due to the relatively large size of the investi- 
gated regions, loop-like intensity enhancements with a size of 
< 1 °were not searched. The original ISSA I^ SA and Ij^ A sur- 
face brightness values were transformed to the COBE/DIRBE 
phot ometric system, usin g the conversion coefficients provided 
bv lWheelock et alJdl994l) : 

- I 60 = 0.87 x lf Q SA + 0.13 MJysr 1 

- Loo = 0.72 x lf^ A - 1.47MJysr- 1 

Dust IR emission maps by ISchlegel et ail Jl998h (SFD) 
were investigated to derive parameters describing our loop fea- 
tures (see Sect. 2.2). The main differences of the SFD 100 /mi 
map compared to the ISSA maps are the following: 

(1) Fourier-destripping was applied, 

(2) asteroids and non-Gaussian noise were removed, 

(3) IRAS and DIRBE 100 /mi maps were combined, preserving 
the DIRBE zero point and calibration, 

(4) stars and galaxies were removed. 

We analysed the radial surface brightness profiles of the loops 
on the SFD 100 /mi map in order to check the effect of the re- 
moval of the sources mentioned above. We also used the SFD 
E(B - V) maps derived from the dust column density maps. 
In the case of these maps the colour temperature was derived 
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Fig. 1. Intensity profile before (top) and after (bottom) back- 
ground removal (gray solid line) with the main derived param- 
eters (I ex , o~ ex , a m , a out , see Sect. 2.2) 

from the DIRBE 100 and 240 /mi maps, and a temperature cor- 
rected map was used to convert the 100 /mi cirrus map to a map 
proportional to dust column density. 

Shape: We approximated the shape of a possible loop by 
an ellipse, which was then fitted using a 2D least-square fit 
method. An ellipse shape is expected from SN or stellar wind 
shells, since (1) non-spherical explosion (wind) may occur, in 
the most extreme case creating a ring, rather than a shell, and 
(2) originally spherical shells are distorted to ellip soidal (i) due 
to the shear in the direction of galactic rotation ( Pa lous et all 
1990) and (ii) d ue to the vertical gravitation al field in the galac- 
tic disc (see e.g. Eh lerova & Palouslll996l) . 

The fitted ellipse is defined with the central (galactic) coor- 
dinates, the minor and major semi-axis of the ellipse, and the 
position angle of the major axis to the circle of galactic latitude 
at the centre of the ellipse. This latter was defined to be '+' 
from East to North (or counter-clockwise). 

Intensity profile: For all of our loops we extracted radially av- 
eraged surface brightness profiles, extended to a distance of 
twice the major (and minor) axis of the fitted ellipse, using 40 
concentric ellipsoidal rings. These surface-brightness profiles 
(ISSA 100 and 60 /mi, SFD 100 /mi and SFD reddening maps) 
were used in the following to determine the basic parameters 
of the FIR emission in the loop. An example is shown in Fig.^ 

2.2. Derived parameters 

Significance: The local background was determined using the 
'non-loop' points in the radial surface brightness profiles, fit- 
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Fig. 2. Examples of loops found in the 2 nd Galactic Quadrant, left: GIRL G134- 28, first described bv lHeilesI i 19841) as GS 137- 
27-17. This supershell is located at a distance at 1.5 kpc and has a size of ~0. 41x0. 3 1 kpc; right: GIRL G l 17+00, supershell 
around the Cas OB5 association in a distance of 2.5 kpc with a size of ~0.3 1 x0.27 kpc ( Moor & KissL 12003). 



ting a 3rd order polynomial. This appropriate background was 
removed from each surface brightness profile points. The in- 
tensity excess I ex was derived as the maximum value of this 
background removed profile. In order to check if this value is 
above the 'noise', we calculated the standard deviation of the 
background removed intensity in the 'non-loop' values, cr ex , 
and defined the significance of the loop as — I ex /cr ex . We 
derived significance parameters on 60 and 100 /im ISSA maps 
and on the SFD lOOjt/m point source removed sky brightness 
and reddening maps. The higher the value of *P is the higher 
the intensity excess of the loop over the background, therefore 
we use this parameter as a 'quality indicator' in the following. 

Relative width: Inner and outer edges of the loop wall along 
the major axis («,-„ and a out , respectively) are defined as the 
radial distance at the full width at half power of the back- 
ground removed intensity profile, I ex /2 (Fig.^ We distinguish 
three regions for a specific loop: (1) loop interior (a<a,„), 
(2) loop wall (a,„ < a < a out ), (3) outer region (a>a OH ,). The 
relative width of the wall of the fitted ellipse is defined as 

W = 1 - din/Clout. 

Colour index: We derived a colour index for our loops 
A/60/A/100, from the radially averaged 60 and 100 fim surface 
brightness profile. This is defined as the slope of the I60 vs. 
I100 scatter plot using the data points of the surface brightness 
profiles in the positions of the loop wall (a,„ < a < a out ) only. 



3. Results 

3.1. The catalogue 

We identified 145 loops in the 2nd Galactic Quadrant. The 
parameters of these loops are summarized in Table A.l. We 




Fig. 3. Distribution of GIRLs in the sky in the 2 nd Galactic 
Quadrant (Aitoff projection), represented by the fitted ellipses 

call these objects 'GiRL'-s, abbreviating 'Galactic JnfraRed 
Loops'. 

The entries of the catalogue are the following: 

1.) Name of the loop, derived from the galactic coordinates 
of the loop centre. The format is : GIRL GLLL±BB, the 
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'GIRL' prefix stands for 'Galactic InfraRed Loops', 'G' 
marks the galactic coordinate system, 'LLL is the galactic 
longitude in degrees, ± is the sign of the galactic latitude 
(+/-), and 'BB' is the absolute value of the galactic latitude 
of the loop centre. 

2. ) Central galactic coordinates of the loop (1 and b) 

3. ) Semi major and minor axis of the fitted ellipse 

4. ) Position angle of the fitted ellipse. The position angle is 

defined as '+' from East to North. Zero position angle is 
pointed to the East in galactic coordinates. 

5. ) Major axes corresponding to the inner and outer edge of 

the loop wall, derived from the 100 /mi surface brightness 
profile 

6. ) ISSA 60 and 100 /mi, SFD 100 //m and reddening signifi- 

cance parameters (*F/6o, T'ioo, 100 and ^ebv, respectively 

7. ) colour index A/60/A/100 of the loop wall 

In the electronic version of the catalogue (URL: 
' http://astro.elte.hu/CFIRLG ') we provide the following 
additional features: 

8. ) 100 /mi ISSA image of the loop with the possible associated 

objects overlaid 

9. ) ISSA 60 and 100 /mi, SFD 100 /mi and SFD E(B-V) surface 

brightness profiles in FITS format 
10.) List of objects apparently associated with the loop 

An associated object has to be placed in the wall or in the 
interior of the loop (defined by concentric ellipses, as described 
above). We considered the following type of possible associ- 
ated objects (references are indicated) 



- dark clouds JDutra & BicaI 1 2002) 

- supernova remnants ( Green, 1 " r ' ' 

- OB-associations (Lan£ 

- pulsars jTavloretallll993l) 

- HII regions (Shar plesslll959i) 

- IRAS point source with molecular core FIR colours 

- IRAS point source with T Tau star-like FIR colours 

Molecular cores and T Tauri stars were selected from the IRAS 
Point Source Catalogue (Joint IRAS Science Worknig Group, 
1988) according to the following criteria: 

- point sources associated with galaxies were excluded 

- photometric qualities are 2 or better at 12, 25 and 60 /mi 

- - molecular cores: 

0.4 < togio(^) < 1.0 & 0.4 < to*io( f£) < 1.3 

- T Tauri stars: 

0.0 < fosio(^) < 0.5 & -0.2 < logmi^j < 0.4 

following the definitions by Emerson (1998). F12, F25 and 
F60 are the 12, 25 and 60 /mi uncorrected IRAS fluxes, re- 
spectively. 

This list of possible associated object does not take into ac- 
count the distances of the individual objects, therefore all ob- 
jects projected to the loop wall or to the interior are included. In 
the case of a few (30) loops, we were able to derive a distance 
for the loop based on the distances of associated objects. This 
is discussed in detail in Sect. 5. The list of associated objects 
is presented in Appendix C. 
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Fig. 4. Distribution of GIRLs in the 2nd Galactic Quadrant, (a) 
Loop-counts vs. galactic longitude. Gray lines with light gray 
filled area: northern galactic hemisphere (b>0°). Dash-dotted 
black lines: southern galactic hemisphere (b<0°). (b) Loop- 
counts vs. the 'sin' of the galactic latitude (equal area belts) 

3.2. Sky distribution and fundamental statistics 

Sky distribution: Fig. |3] presents the distribution of the far- 
infrared loops found in the 2nd Galactic Quadrant. The dis- 
tribution of GIRLs seems to be rather complex. However, there 
are slightly more loops at low (b < 30°) galactic latitudes than 
at high ones (b>30°). 

It is obvious from Fig|4^ that the two galactic hemispheres 
represent different loop distributions. We have tested the hy- 
pothesis that the galactic latitudes of the loop centres on the 
Northern and Southern hemispheres are drawn from the same 
distribution. A Kolmogorov-Smirnov test showed that the ho- 
mogeneous distribution hypothesis could not be accepted at the 
80% confidence level. A test performed for the galactic longi- 
tude distribution on the two hemispheres (Fig|4j)) yielded the 
same result. 

There is a hint in Fig. |5J that most of the loops are part 
of a structure with scales significantly larger than the individ- 
ual loop diameters. The most remarkable feature of this kind 
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is a huge chain of loops seen in the 1»;90°-130°, b=-35°- 
+5°region. In the ranges 1 * 130°-140°, b * +45°- +65°nine 
loops can be found squeezed to a relatively small area, which 
correlates well with the position of the Ursa Maior molecu- 
lar complex connected to th e expanding shell of the NCP loop 
JPound A Goodman! fl997lh 



40 F 




r e ff (cleg) 

Fig. 5. Apparent size distribution of GIRLs. Size is presented in 
effective radius, r e ff, the geometrical mean of semi-minor and 
-major axes of the fitted ellipse 



Size distribution: We defined the effective radius of our loops 
as the geometrical mean of the semi-minor and -major axis of 
the fitted ellipse, r e ff = Vfl X b. For the whole sample of loops 
we obtained m(r e ff) = 3°2 and E(r e ff) - 2°7 for the sample 
mean and expectation value of the effective radius, respectively. 

The size distribution is presented in Fig.|5] Majority (78%) 
of the loops are in the apparent effective radius r e ff <4°, and 
there are only a few (9) loops r e yy > 7°. The lower end of the 
size distribution is limited by our detection method, which ex- 
cludes features with a size smaller than ~l°in diameter. 

Distribution of relative width: Fig. |6] shows a Gaussian- 
like distribution, with an average value of (W) = 0.29 and a 
FWHM = 0.04 for the relative width W of the loop walls. 

Axis ratio distribution: As mentioned in Sect. 1, it is widely ac- 
cepted, that galactic loops - identified in any wavelength range 
- are 2D projections of 3D cavities in the ISM. We compared 
the expected distribution of the ratio of minor to major axes of 
the sky-projection of 2D ring to the distribution of the observed 
minor to major axes ratios in Fig.Q 2D rings are assumed to be 
randomly oriented in <p and 9 polar angles. This figure clearly 
shows a significant difference between the two distributions. It 
is not likely that most of our objects are just the projections of 
2D rings. 
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Fig. 6. Distribution of the width of the loop walls relative to the 
effective radius of the loop (W = 1 - a,„/a „,) 
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Fig. 7. Comparison of the expected distribution of minor to 
major axis ratio (assuming randomly oriented 2D rings, solid 
line) with the distribution of the observed axis ratios (gray his- 
togram). 

Significances: In Fig.|9]we compare the significances derived 
from ISSA 100 /mi, ISSA 60 /an, SFD 100 //m and SFD red- 
dening maps, ¥/ioo, f/60, ^Sfdiqo, ^ebv (see also Table A.l). 
The Pearson correlation coefficients are listed in Table [0 The 
^noo, ^m, ^sioo and *¥sebv significances are well correlated, 
and in all 145 cases resulted in a definite confirmation of the 
suspected ellipsoidal structure. 

Fig. |8] presents the distribution of ISSA 100 /an and SFD 
E(B-V) significances (T^ioo and ^ebv, respectively). Based on 
the definition of *P (see Sect 2.2) one can set a limit of *P = 3 as 
a value separating weak and prominent loops (prominent loops 
are above the "3cr" level of background fluctuations). However, 
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correlated 


correlation 


significances 


coefficients 


™J 100 - f/60 


0.77 


"P/ioo - IS 100 


0.82 


^/lOO - ^SEBV 


0.68 


100 - ^SEBV 


0.79 



Table 1. Linear Pearson correlation coefficients of signifi- 
cances derived from ISSA 100 yum OF/ion), ISSA 60 fjxa (T /60 ), 
SFD 100 //m OFsioo) and SFD E(B-V) Q¥ S ebv) maps. 
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Fig. 8. Distribution of the "F/ioo (gray solid line) and ^ebv 
(black dash-dotted line) significance values 



as seen in Fig. [8] the majority of our loops show high sigma 
values, regardless of the significancy type OF/ioo, *F/60, "Fsioo 
or ^ebv) an d on ly ver Y f ew loops have *F < 3. 

Bright fore- or background sources not related to the loop 
may cause low significances. An example is GIRL G125+09, 
which has low *Fnoo and *Psioo (2.4 and 2.6, respectively) but it 
is outstanding on the reddening maps C^ebv = 6.0) and is very 
conspicuous at 60 fim (f/6o = 17.2). I n addition, this lo op is 
clearly detected in the HI 21 cm line dKiss et all \200(j . The 
main reason behind the low significances at 100 fim is the pres- 
ence of an extremely bright region in the south-eastern part of 
the loop related to NGC 7822, S 171 and to the CepOB4 asso- 
ciation, which dominates the background of GIRL G125+09. 



Far-infrared colour of the loop walls: As seen in Fig.fTUlmost 

of the loops have a far-infrared colour similar to that of the 
extended galactic background emission or galactic cirrus with 
an average value of Aleo/AIjoo = 0.25+0.12. In the case of a 
few loops the far-infrared colour significantly deviates from 




10 
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103 



Fig. 9. Comparison of significances extracted from maps of dif- 
ferent data sets, (a) ISSA 100 /an vs. ISSA 60 /mi; (b) ISSA 
100//m vs. SFD 100 //m; (c) ISSA 100 yt/m vs. SFD reddening 
maps; (d) SFD 100 fim vs. reddening 

that of the extended background, especially in the vicinity of 
the galactic plane (|b| < 10°, see Fig. I10> . where loops are 
"warmer" according to their FIR colours. A rather high scat- 
ter of FIR colours with high error bars can be observed for the 
highest galactic latitudes, which could also be related to the 
less reliable accuracy of the surface brightness calibration in 
these fields. We conclude that the loops describe the structure 
of cirrus in the few to ~20° scales. 

Excess extinction in the loop walls: Using the radially av- 
eraged SFD E(B-V) (reddening) profiles extracted from the 
appropriate maps we estimated the excess visual extinction 
of the loop walls. Visual extinction values were calculated as 
Ay = R x E(B - V), where R is the ratio of total over selective 
extinction and E(B-V) is the reddening value in a specific sky 
position. R was choosen to be uni form and set t o the standard 
value of 3. 1 of the diffuse medium riKriigell 12003). The average 
value of Ay was derived for the inner-, wall- and outer regions 
(see Sect. 2.2), separately for all loops (A™, A™" and A°"', re- 
spectively). We also derived an excess visual extinction value 
for the loop wall, which was defined as 



AAv 



2- A 



wall 



A^ + A""' 



(l) 



AAv values show a clear dependence (despite the relatively 
large scatter around the main trend) on the galactic latitude 
(Fig. II lh). width values <07 1 01 at high galactic latitudes, while 
reaching AAy » 10'" in the vicinity of the galactic plane. 

The relative strength of the excess extinction AAy relative 
to the background extinction (A^ 8 = | [A" + A° ut ]) is usually 
low (in average ~10%, as shown in Fig.fTTb) except in the case 
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Fig. 10. Distribution of the colour indices of the loop walls (Aleo/AIioo) vs. galactic latitude. Error bars in galactic latitude 
represent the extension (i e ff) of the loop. Error bars of the colour indices are the formal errors of the Iioo vs.l6o scatter plots. The 
thick gray curve represents the main trend. 



of a few loops near the galactic plane. This corresponds to a 
~10% increase of column density in the loop walls. 

3.3. Distances to individual loops 

Derivation of the distance: Many of the major physical prop- 
erties of a loop (mass, size, particle number density, etc.) can 
only be tackled if the distance is known. In order to estimate the 
distance of the loops, we used associated objects with known 
distances. There are two major classes of these objects: (1) 
They may be responsible for the formation of the loop itself, 
or are the remnants of the loop-creating events (we considered 
pulsars, luminous stars and/or their associations; see Sect. 13. II 
for references). (2) They are part of the loop, but identified in 
anot her wavel ength range (dark clouds and molecular clouds, 
butra & BicaLl2002h . 

We applied the following main criteria when selecting dis- 
tance indicators: 

i) Luminous stars had to have the power to form the loop, 
therefore only O stars and B stars with earlier subtype (or 



their associations) were taken into account. The sky posi- 
tion of the stars must be well confined in the central part of 
the loop (although in some cases fast proper motions might 
have placed the star far from its original position). 

ii) There must be at least two objects (luminous star, pulsar, 
dark/molecular cloud) with similar known distances associ- 
ated with the loop, otherwise the distance hypothesis would 
be rejected. 

iii) The distances should fit each other within their uncertain- 
ties. 

Galactic distribution: Based on the distances derived for some 
of our loops we drew the 3-dimensional distribution of GIRLs 
with known distances in Fig. 1121 

Since our distance indicators are mostly objects which are 
related to the spiral arms of the Galaxy, the distribution reflects 
the position of the spiral arms, introducing a selection effect. 
Loops closer than 800 pc are related to OB associations and 
luminous stars of the Local Arm and to the nearby interarm 
region, while another well-defined group can be identified at a 
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Fig. 11. Excess visual extinction in the loop walls derived from 
the SFD E(B-V) surface brightness profiles, (a) galactic lat- 
itude vs. excess visual extinction AAy (b) galactic latitude 
vs. the ratio of the excess visual extinction AAy to the av- 
erage of the inner and outer visual extinction values (A^ 
i [A » + A° ut ]) 
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distance of 1.8-2.7 kpc, in agreement with the distance of the 

J 1 1 1 

Perseus Arm in that direction ( HumnhrevSL 119781) . There is a 

hint that the majority of the GIRLs are within a distance of 

~800pc. The most prominent loops remained observable up to 

a distance of ~3 kpc. 




(po) 




Fig. 12. Loops with known distances in the 2nd Galactic 
Quadrant, projected to the galactic plane, in cylindriclateal co- 
ordinate system. The top panel shows the loops with projected 
distances up to 4 kpc. The bottom panel is a magnified view of 
the region closer than 1 kpc. Different gray tones in both the 
upper and lower panels code the vertical distance of the loop 
centres from the galactic plane (see colour bar). 



4. Discussion 

Size distribution: In the Galaxy most of the ISM is con- 
fine d to a thin disc, formin g the Galactic Disc itself (see e.g. 
Dickev & Lockma il l 199(1 for a review). Therefore it is ex- 
pected that any object formed in the ISM should somehow 
reflect the distribution of the originating medium. Accepting 
that our loops are features of the cold diffuse ISM, these have 
been formed in the galactic disc with a scale height of ~250pc. 
Therefore loops at high galactic latitudes are close to our Solar 
System. On the other hand, in an exponential disc the num- 
ber density of the gas and dust is decreasing rapidly with the 
increasing height from the galactic plane. This allows the gen- 



eration of larger SN- or stellar wind bubbles even possessing 
the same energy input as the ones closer to the galactic plane. 
Thus one would expect to observe loops with a large appar- 
ent size at high galactic latitudes. Still, we found no correlation 
between the galactic latitude of the loop centres and the effec- 
tive radii, neither on the northern, nor on the southern galactic 
hemisphere. 

The sky distribution presented in Figs. [3] and |4] shows 
clearly that the distribution of the GIRLs found at the 2 nd 
Galactic Quadrant is far from that expected e.g. from the expo- 
nential disc hypothesis, in which the regions near to the galactic 
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Fig. 13. (a) Distribution of distances among the loops presented 
in Fig^](solid lines with gray filled area). Dashed-dotted lines 
represent the distribution of distances projected to the galac- 
tic plane, (b) Distribution of effective diameters (2R e ff) in the 
same sample, (c) The same as b but displayed and binned on a 
logarithmic scale. A relation of log2R e fj-logN was fitted for 
1 .4 < 2R e /f < 2.6. The points used for the fit are marked by as- 
terisks (see Sect. 4 for details). 



plane should have the largest number of loops and this number 
should decrease drastically by the increasing galactic latitude. 
The relatively low number of loops in the galactic midplane 



may be explained by multiple overlapping of the loop walls 
projected on the sky, i.e. many loops might be undiscovered 
due to confusion. Intervening molecular clouds may also con- 
tribute to this confusion. 

As discussed in Sect. 3.3, there are two separate distance 
groups in the known-distance sample. Accepting that the ma- 
jority of our loop are located in the nearby group (discussed in 
detail in the next paragraph) at a distance d< 800 pc, one can 
obtain the expected ratio of loop counts above and below \b\ - 
30° for a random distribution with a scale height of zh = 250 pc 
(Nakanishi & Sofuel 120031) . This results in (Ni ow /Nhi g h) ca ic = 
2.15 which should be compared with the observed ratio of 
{Niow/Nhi g h) bs = 1 .28 (without the contribution of the Perseus- 
Arm group). If this discrepancy is caused by confusion as dis- 
cussed above, we may miss ~40% of the total count in the 
galactic midplane. 

It should be noted that molecular gas extends to galactic lat- 
itudes higher (both at the northern M artmann et all 1 19981 fig. 
2) and the southern ( Magnan i et alll2000l fig. 1) galactic hemi- 
sphere) in the 2 nd Galactic Quadrant than in the other parts of 
the sky. This anomalous feature may also explain the richness 
of high galactic latitude loops in our sample. 

The large number of loops at high galactic latitudes might 
be a selection effect (loops are easier to be identified due to the 
lack of confusion). Still, one needs a process that can generate 
loops far from the galactic plane. One of these processes can 
be the infall of clouds from the halo, and their collision with 
the galactic disc. We kno w some examples, e.g. the NCP-loop 
(Me verdierks et allll99lh . but the large number of such loops 
would require a high rate of cloud infall, which is not supported 
by observations jEhlerova & PalousU 1996). 

Another possibility is that turbul ent processes in the ISM 
create ' loop-like ' structu res JKlessen et all Eoool iKorpi et all 
1999; Wada et al., 2000). Nonlinear development of turbulent 
instabilities may form voids in the ISM without the need of 
stellar energy injection. These instabilities may be responsible 
for the large number of loop features at high galactic latitudes 
since these parts of the sky lack suitable stellar sources (mas- 
sive O and/or B stars). From the observational point of view the 
'loops' created by these two different processes are practically 
indistinguishable, and therefore we are not able to determine 
the relative importance of these processes. 

It is also clear from the figures that local effects play a very 
important role forming the distribution. We can mention e.g. 
the high latitude group at 1= 130°-140°, b=45°-65°, which 
is related to the huge NCP-loop and also belongs to a group of 
high latitude molecular clouds. Similar groups contribute much 
to the peaks in Fig.[4k and b, and might be responsible for the 
deviations from the expected distributions. Since these struc- 
tures are relatively nearby (e.g. the distance of the NCP-loop is 
~330pc), even if we can detect loops up to a distance of ~3 kpc, 
the main contribution comes from loops at small distances (see 
Sect. 3.3). 

In Fig.ll3t we displayed the size distribution on a logarith- 
mic scale and fitted the large diameter part (data points marked 
by asterisks) with a power-law (linear in the log-log plot). 
Simila r investigation has been carried out by lOev & Clarke! 
i 19971) for the Small Magellanic Cloud, Holmberg II, M31, 
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M 33 and by iKim et all d2003l) for the HI data of the Large 
Magellanic Cloud. We found a power law index of s = -0.47 ± 
0.11 (following |Kirnetal.Ll2003l) or^n = 1 .24 + 0.30 using the 
transformation bv lOev & Clarke (1997), which is lower than 
that of the other galaxies mentioned above. We note that the 
distribution has a large scatter and it is better described by a 
double-peak one than by a power-law. 
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Fig. 14. Relationship between distance and effective diameter 
(2R e /f) for loops with known distances (see Table B.l). The 
shaded areas represent the apparent size limits 2r e /f < 1° and 
d>40°. Due to the limitations discussed in Sect. 2.1 we were 
not able to detect far-infrared features of these sizes. 

Correlation of distance and size: A good correlation can be 
found between the physical size of the bubble and the derived 
distance (for the sample in Sect. 3.3) when we derive linear 
Pearson correlation coefficients for the whole available dis- 
tance range, up to 3.5 kpc (0.66 and 0.73 for the linear-linear 
and log-log plot, respectively). Fig. ^] clearly illustrates the 
limitations of our study: due to the limitations in apparent size 
(l°<d<40°) we identified loops in the allowed region only. 
This selection effect is more severe for the small-sized part of 
the distribution (1° <2r e f/), while even the largest loop is a 
factor of ~2 smaller than the largest detectable feature. The 
distribution shows two groups, one at D <800pc and one at 
D >800pc (Fig.ll3h). Within the nearby group practically there 
is no correlation between the distance and physical size (with a 
linear correlation coefficient of 0. 17), therefore this group is not 
affected by the selection effect. Accordingly, we may derive a 
statistical result for the local ISM, while the Perseus Arm group 
is strongly biased. The size distribution of the nearby group 
lacks the large bubbles (2R e tt > lOOpc) which are observable 
in the Perseus-group. The lack of large bubbles in the Local 
Arm group might be the consequence of a non-continuous star- 



formation history in the Milky Way or related to the finite ex- 
tension of the Local Arm. 

Volume filling factor: The presence of the Perseus arm group 
is probably restricted to the galactic plane and represent only 
a small subsample. Therefore it is plausible to assume that 
the majority of the loops with unknown distances belong to 
the nearby (D < 800 pc) group. Applying the distance - size 
distribution obtained for the nearby group to all loops with 
unknown distances (except the loops of the Perseus group), 
it is possible to give a rough estimate for the volume filling 
factor (or porosity Q, ICox & Smith! 1 19741) of our loop interi- 
ors in the galactic disc. Assuming that the loops/bubbles were 
formed in the cold neutral interstellar medium (justified by 
their FIR colours) using zh = 250 pc for the scale height of HI 
dNakanishi & Sofud.120031) in the solar neighbourhood, we ob- 
tain a volume filling factor of f 2 nd ~ 4.6% for the 2 nd Galactic 
Quadrant for D<800pc. Taking into account an additional 
~40% of loops in the midplane (see above) an upper limit of 
f 2 nd x>6.4 % can be obtaine d. T hese values confirm the filling 
factors bv lFerrierel i 19981) and iGazol-Patino & Passotl fl999) 
(/ = 3-8% for 8.5 <R e < lOkpc galactocentric distances). 

5. Summary and outlook 

We summarize the main outcomes of our paper below: 

- One fourth of the far-infrared sky was surveyed and 145 
loops were identified. 

- The loop sample is representative in the Local Arm 

- The FIR loops outside the midplane possess ~10% column 
density excess and a cirrus-like A/60/A/100 ~ 0.25 ratio 

- The midplane loops hold several times denser walls and a 
far-infrared colour of Afe)/ A/100 > 0.4 

- Distance estimates are given for 30 loops which trace the 
cold ISM in the Local Arm and in the Perseus Arm 

- We estimated the hot gas volume filling factor to be 
4.6% <<f 2 *d <6.4% in the Local Arm in the 2 nd Galactic 
Quadrant 

As a continuation of this study the 1 st , 3 rd and 4 th galactic 
quadrants are surveyed as well. The complete sample and a 
comparison with HI distribution will be given in a forthcoming 
paper. 
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